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The stellar populations of galaxies hold vital clues to their formation histories. In this paper we present results 
based on modeling stacked spectra of early-type galaxies drawn from the Sloan Digital Sky Survey (SDSS) as a 
function of velocity dispersion, a, from 90 km s _I to 300 km s -1 . The spectra are of extremely high quality, with 
typical signal-to-noise ratio, S/N, of 1000 A -1 , and a wavelength coverage of 4000A-8800A. Our population 
synthesis model includes variation in 16 elements from C to Ba, a two-component star formation history, the 
shift in effective temperature, A T e ff , of the stars with respect to a solar metallicity isochrone, the stellar initial 
mass function (IMF), amongst other parameters. In our approach we fit the full optical spectra rather than a 
select number of spectral indices and are able to, for the first time, measure the abundances of the elements 
V, Cr, Mn, Co, and Ni from the integrated light of distant galaxies. Our main results are as follows: 1) light- 
weighted stellar ages range from 6-12 Gyr from low to high cr; 2) [Fe/H] varies by less than 0. 1 dex across the 
entire sample; 3) Mg closely tracks O, and both increase from w 0.0 at low a to ~ 0.25 at high cr; Si and Ti show 
a shallower rise with cr, and Ca tracks Fe rather than O; 4) the iron peak elements V, Cr, Mn, and Ni track Fe, 
while Co tracks O, suggesting that Co forms primarily in massive stars; 5) C and N track O over the full sample 
and [C/Fe] and [N/Fe] exceed 0.2 at high a ; and 6) the variation in A T e ff with total metallicity closely follows 
theoretical predictions based on stellar evolution theory. This last result is significant because it implies that we 
are robustly solving not only for the detailed abundance patterns but also the detailed temperature distributions 
(i.e., isochrones) of the stars in these galaxies. A variety of tests reveal that the systematic uncertainties in our 
measurements are probably 0.05 dex or less. Our derived [Mg/Fe] and [O/Fe] abundance ratios are 0.05-0.1 
dex lower than most previous determinations. Under the conventional interpretation that the variation in these 
ratios is due to star formation timescale variations, our results suggest longer star formation timescales for 
massive early-type galaxies than previous studies. Detailed chemical evolution models are necessary in order 
to translate the abundance ratio distributions of these galaxies into constraints on their formation histories. 
Alternatively, these data may provide useful constraints on the nucleosynthetic pathways for elements whose 
production is not well understood. 

Subject headings: galaxies: stellar content — galaxies: abundances — galaxies: early-type 



1. INTRODUCTION 

We now have a working theory for the evolution of large 
scale structure in the universe thanks to the largely successful 
paradigm of a universe dominated by cold dark matter. In this 
paradigm cosmological structures form "bottom-up" — small 
structures collapse early and the massive clusters assemble 
late. One of the more interesting puzzles in this paradigm is 
the realization that the most massive galaxies in our universe 
seem to have formed their stars earlies t, with lower-mass 
galaxies forming t heir stars later (e.g., TrageretaTJ 119981; 
Thomas et al. 2005). This phenomenon goes by many names, 
including "downsizing" and "anti-hierarchical growth". No 
ab initio model for galaxy formation is able to reproduce these 
basic features of the observed galaxy population. The appar- 
ent tension between the bottom-up nature of cold dark matter 
and the top-down star formation histories of galaxies can be 
reconciled by recognizing that the formation time of the stars 
need not be related to the assembly time of the stars into their 
present day galaxies. 

From an empirical point of view, the formation and evolu- 
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tion of galaxies can be probed via two general techniques. The 
first is through lookback studies where one observes, statisti- 
cally, the progenitors of present day galaxies at progressively 
higher redshifts. The second is through studying the present 
day properties of galaxies, including their stellar populations, 
structure, and kinematics, in order to learn about their past 
evolution. This latter technique is often referred to as the 
archeological approach, and is the subject of this paper. 

It has long been appreciated that optical through near- 
infrared (NIR) spectra of old stellar systems are rich in ab- 
sorption features that provide clues to the ages, metal licities, 
and a bundance patterns of the stars (for a review, see IConrovl 
120131) . Since the 1980s the standard analysis technique has 
been to measure and model the strength of ste llar absorption 
features through the Lick /IDS index system (Burstei n et al.l 
U984tlWorthey et al.lli994l) . Authors typically model a hand- 
full of these indices in order to estimate various stellar popu- 
lation properties. 

Analysis of the Lick indices has revealed that early-type 
galaxies are enhanced in the a element Mg compar ed to the 
abund ance patterns of stars in the Galactic disk (IWortheyl 
1994). The [a/Fe] ratio has long been known to be sensitive 
to several parameters, including the timescale of star forma- 
tion, the slope of the initial mass function (IMF) at > 1M , 
the delay time distribution of Type la superno vae (SNe), and 
the preferential loss of metals via winds (e.g., Tinsley Hl979l : 
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Figure 1. Model spectrum for an age of 13 Gyr and solar metallicity. The spectrum has been smoothed with a velocity dispersion of a = 350 kms~' , equal to the 
smoothing applied to the early-type galaxy data analyzed in this paper. Strong features are labeled. Also included is the location of the true stellar continuum, 
which is the spectrum that would be observed in the absence of all line opacity. In this figure the model spectrum is computed entirely from synthetic stellar 
spectra, whereas for the main analysis the synthetic spectra are only used differentially. 



iThomas et alj [l999). The sensitivity of [a/Fe] to all of these 
processes significantly complicates the interpretation of this 
ratio. Nonetheless, it has become standard practice to as- 
sume that the [a/Fe] ratio is telling us something about the 
timescale of star formation. By c omparing to simple c losed- 
box chemical evolution models, Thomas et al. (2005)) con- 
cluded that the [a/Fe] ratios of the most massive galaxies 
in their sample, with a w 300 kms -1 , implied star formation 
timescales of only ~ 200 Myr. These short inferred star for- 
mation timescales, if correct, would have profound implica- 
tions for our understanding of the formation of the most mas- 
sive galaxies in the universe. 

The analysis of spectral indices is beginning to give way 
to the modeling of the full optical-NIR spectra of stel- 
lar systems. This approach has great potential, and is 
now frequently used to measure st ar formation histories and 
metallicities non-parametrically dCid Fernandes et all 120051; 
lOcvirk etal .1120061; iToieiro et al.l 120091) . Recently, such full 
spectrum fitting models have been exten ded to include varia- 
tion in the elemental abundance patterns dWalcher et al.l2009t 
iConrov & van D okkum 2012a), paving the way to measure 
many more parameters from integrated light spectra than 
would ever be possible with spectral indicesj. 

This paper is part of an ongoing series focused on modeling 
very high quality optical-NIR spectra of early-type galaxies 
with full-spectr um fitting stellar population sy nthesis (SPS) 
techniques . In IConrov & van Dokkuml (I2012al CvD 12), we 
presented the model, and in IConrov & van Dokkuml (12012b!) 
we used that model to measure the low-mass IMF from the 
spectra of 34 nearby early-type galaxies and the nuclear bulge 

4 Indices typically have a central bandpass of roughly 30 A wide, and 
so measuring indices from spectra is roughly equivalent to heavily boxcar 
smoothing the data. In addition, indices are sensitive not only to the fea- 
ture^) within the central bandpass, but also the features in the pseudocontin- 
uum, and this can significantly complicate interpretation of the indices. Full 
spectrum fitting does not suffer from these limitations. 



of M 3 1 (b ased on data present ed in Ivan Dokkum & Conrovl 
120121) . In IConrov et al.l (120131) we presented results on the 
neutron-capture elements Sr and Ba based on our local galaxy 
sample and on the stacked spectra that are the focus of this 
work. In the present paper we measure the ages, detailed 
abundance patterns, and effective temperatures of the stars in 
early-type galaxies drawn from the SDSS. In a subsequent 
paper we will discuss the low-mass IMFs inferred for these 
galaxies. 

The rest of this paper is organized as follows. Section|2]pro- 
vides an overview of the SPS model and our fitting technique, 
Section [3] describes the data, and Section |4]presents a test of 
the model by fitting to spectra of metal-rich star clusters. Our 
main results are presented in Section and in Section [6] we 
compare our derived properties to results based on other mod- 
eling techniques. A discussion and summary of our results is 
provided in Sections [7] and [8] 

2. MODEL & FITTING TECHNIQUE 

The SPS model u sed he rein was developed in 
IConrov & van Dokkuml d2012al CvD 12). The model 
adopts standard SPS techniques, including constructing 
libraries of isochrones and stellar spectra. Empirical spectra 
form the core of the model. We make use of two separate 
libraries, the MILES library, which covers the wavelength 
range 0.35fj,m < A < 0.74fim (ISanchez-Blazquez et alJ 
2006), and the IRTF library of cool sta rs, which covers the 
wavelength range .8 \jim < A < 2.4/im (Cushin g"et al.ll2005l 
iRavner et alj|2009l) . The empirical stars are of approximately 
solar metallicity and have solar abundance patterns. We 
have computed a large grid of model stellar atmospheres 
and spectra in order to construct response functions (i.e., the 
relative change in the spectrum of a star due to a change in 
the abundance of a single element). We chose 20 positions 
along a 13 Gyr isochrone to compute synthetic spectra for 
each abundance pattern. The model atmospheres and spectra 
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Figure 2. Relative response of the blue spectral region to an increase in the abundance of nine different elements. Abundances have been increased by 0.3 dex 
for all elements except for C, which is increased by 0.15 dex. The changes are with respect to a solar metallicity model with a MW IMF and an age of 13 Gyr. 
The models have been smoothed to a velocity dispersion of a = 150 km s - '. Notice the different y-axis range in each panel. It is evident from both this figure and 
Figures |3]and|4]that essentially every spectral region is influenced by one or more elements. It is therefore very difficult to define a spectral index that is sensitive 
to only one element. 



were computed with the AT LAS model a t mosph ere and 
spectrum sy nthesis package (iKuruczl 119701 11993b . ported 
to Linux by ISbordone et al.l (120041) . Specifically, we use 
the ATLAS 12 code, computing new atmospheres for each 
change in abundance pattern. The line list was provided by R. 
Kuruc20, including linelists for TiO and H2O, amongst other 
molecules. The spectral response functions were applied 
to the empirical stellar spectra in order to create models 
with arbitrary abundance patterns. Notice that the synthetic 
spectra are only used differentially. The model allows for 
arbitrary variation in the IMF and spans ages from 3-13.5 
Gyr. See CvD12 for further details regarding the model. 
Figure Q] shows the model spectrum for a 13 Gyr, solar 
metallicity population, smoothed to a velocity dispersion of 
cr = 350 kms -1 . 

There have been several updates to the CvD12 model. Most 
importantly, the synthetic atmospheres and spectra have been 
recomputed with an updated line list kindly provided by R. 

5 kurucz . harvard . edu/ 



Kurucz. This line list includes many more iron lines and 
updated wavelengths and hyperfine splitting of many lines. 
The model now contains variation in many more elements, 
described below. We have also supplemented the MILES 
empirical library with stacked SDSS M dwarf spectra from 
iBochanski et al.l (120071) . This turns out to be important be- 
cause the MILES library contains only a few M dwarf spec- 
tra, and none cooler than M6V. The SDSS M dwarf template 
spectra cover the entire M dwarf spectral sequence. The fi- 
nal update since CvD12 is the inclusion of emission lines 
in the fitting (rather than masking regions of potential emis- 
sion line contamination, as done previously). Currently we 
include 13 emission lines, including individual components 
of [O n]A3727, [Ne m]A3868, US, H 7 , H/3, [O m]A5007, 
[N i]A5201, [N n]A6547,6583, Ha, and [S 11] A67 16,6731. 
We include separate parameters for the strengths of each of 
these lines (i.e., we do not enforce particular line ratios, even 
amongst doublets with known, constant line ratios). The in- 
trinsic line shape is assumed to be a gaussian with a single 
width for all lines. This width is also included in the fit. 
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Figure 3. Same as Figure[3] now showing the red spectral region. Abundances 
by0.15dex 

The effect of varying elemental abundances on the model 
spectra are shown in Figures |2] [3] and |4] These figures show 
the response of a 13 Gyr, solar metallicity model to an in- 
crease in the abundance of a single element. We emphasize 
that the model atmospheres and spectra are always treated 
self-consistently (i.e., the elemental concentration is increased 
both in the computation of the model atmosphere and of the 
model spectra). In many cases the effect on the spectrum 
from a change in an element is straightforward, e.g., increas- 
ing [Mg/H] results in stronger Mg I lines, but other effects are 
more subtle. For example, because Mg is a major electron 
donor, increasing [Mg/H] affects the ionization equilibrium 
and therefore affects many species of other elements. For ex- 
ample, it is for this reason that [Mg/H] affects the strength of 
the Ca II lines at ~ 0.86/xm. 

We follow iConrov & van Dokkuml (1201 2bh in fitting the 
model to data. In its present form the model contains 40 free 
parameters, including the redshift and velocity dispersion, a 
two-part power-law IMF, two population ages, four nuisance 
parameters, 13 emission line strengths, the velocity broaden- 
ing of the emission lines, and the abundances of C, N, Na, 
Mg, Si, Ca, Ti, V, Cr, Mn, Fe, Co, Ni, Sr, and Ba, and O, 
Ne, S (the latter three are varied in lock-step). These parame- 



have been increased by 0.3 dex for all elements except for C, which is increased 

ters are fit to the data via a Markov Chain Monte Carlo fitting 
technique. The data and models are split into four wavelength 
intervals (described below) and, within each interval the spec- 
tra are normalized by a high-order polynomial (with degree 
n wher e n = (A max - A m j n )/100A; see lConrov & van Dokkuml 
2012b, for details). For the purposes of the present article, 
the emission line strengths are treated as additional nuisance 
parameters, and in many cases their strengths are constrained 
to be very low, which is not surprising given the sample defi- 
nition described in the next section. 

One of the nuisance parameters included in the model is a 
shift in effective temperature, A T e ff , applied to all stars used 
in the synthesis. As described in CvD12, this parameter is 
meant to capture changes in the isochrones due to changes in 
abundance patterns of the model stars. The motivation for 
this shift is shown in Figure [5] In the top panel we show 
isochrones obtai ned from the Dar tmouth Stellar Evolution 
Database (DSEP; lDotter et al.f2 008) for a range in [Fe/H] and 
[a/Fe] . In the bottom panel we have shifted each isochrone by 
a constant amount in T e a until a match was achieved with the 
solar metallicity case. Notice that in all cases a constant shift 
does an excellent job of accounting for modest variation in 
abundance patterns from the solar metallicity model. More- 
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Figure 4. Same as Figure[2] now showing the elements encompassing the iron peak (V, Cr, Mn, Fe, Co, Ni). The change in the model spectra have been computed 
for a 0.3 dex increase in the abundance of each element. 



over, it is interesting that the isochrone with [Fe/H]= -0.1 
and [a/Fe]= 0.2 is almost identical to the solar metallicity 
isochrone. Thus, over the metallicity range of interest in this 
work, a simple shift in T e f[ is sufficient to account for el- 
emental abundance effects on the isochrones. Because the 
abundances of many elements are varied simultaneously in 
the model, we use the free parameter AT e ff to encompass 
their combined effect, rather than attempting to compute self- 
consistent isochrones from first principles. 

3. DATA 

3.1. SDSSData 

The main analysis presented here results from fitting our 
updated CvD stellar population model to very high-S/N 
stacks of passive galaxies, binned by stellar velo city disper- 
sion. These galaxies are selec ted from the SPSS (lYork et al.l 
120001) Main Galaxy Survey fStrauss et al.1 12002b Data Re- 
lease 7 (lAbazaiian et al.l 12009b . within a narrow redshift in- 
terval (0.02 < z < 0. 06). Following the methodology of 
Peek & Gravesl d20l0h . we select passive galaxies by re- 
quiring that they have no detected emission in Ha nor in 
[O n]A3727. 

The individual spectra have moderate S/N, typically ~ 20 
A" 1 , so we stack the spectra of hundreds of galaxies to achieve 
the very high S/N needed for detailed abundance analysis. 
Passive galaxy star formation histories vary strongly as a func- 
tion of their stellar velocity dispersion, but also at fixed a de- 
pen ding on how f ar the y scatter off th e Fundamental Plane 
(FP: iGraves etaniMTlTspringob et al.ll20"l2h . We therefore 
have chosen to stack galaxies in bins of a, in order to study 
the main dimension of stellar population variation, but restrict 
our sample to galax ies that he on the FP (i.e., in the central 
FP slice defined in lGraves et al.ll2()10l) . The two-dimensional 
space of galaxy star formation histories will be explored in a 
future paper. 



We stack the spectra in seven bins of velocity dispersion, 
with mean values of a = 88, 112, 138, 167, 203, 246, and 
300 km s -1 . The effective resolution of the data therefore 
varies from R k, 400- 1200 owing to intrinsic Doppler broad- 
ening. The SDSS spectra are obtained through 3" diameter 
fibers, which sample the inner ps 0.8 R e for the smallest a bin 
and the inner ps QAR e for the largest a bin. Each individ- 
ual spectrum was continuum-normalized and convolved to an 
effective dispersion of 350 km s" 1 before stacking, and each 
spectrum contributed equally to the stack. Problematic pixels 
(e.g., those under bright sky lines) were masked before stack- 
ing. The resulting stacked spectra have S/N at 5000 A that 
ranges from ps 500 A" 1 to ps 1800 A" 1 . Four wavelength in- 
tervals were defined for the SDSS stacks, 0.40/^m-0.48/.tm, 
0.48/im-0.58/im, 0.58^m-0.64^m, and 0.80/im-0.88^m for 
fitting models to data. 

In Figure [6] we show the stacked spectra in three a bins. 
The fluxes have been continuum normalized in order to focus 
attention on the differences in the line-strengths between dif- 
ferent stacks. Notice that the differences are relatively subtle, 
with features strengths varying by only a few percent from the 
lowest to the highest dispersion bins. 

3.2. Star Cluster Data 

In this paper we will also analyze integrated light spectra 
of metal-rich star clusters in order to test the model. We fo- 
cus on two ~ Z Q clusters, M67 and NGC 6528. The former 
is an intermediate-age open cluster with solar-scaled abun- 
dance ratios, while the latter is an old bulge globular cluster 
with enhanced abundance ratios. For M67 we adopt an age 
of 4 Gyr derived from its color-magnitud e diagram (CMD; 
iVandenBerg e"t1II 120071 iMaeic et alJl20Toh . [Fe/H]=0.0, and 
solar-scaled abundance ratios with an uncertainty of . 1 dex 
dTautvaisiene et al.ll2000t IShetrone & Sandquisj 12000). For 
NGC 6528 we adopt an age of 11 Gyr, consistent with re- 
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Table 1 

Results from Modeling Stacked SDSS Early-Type Galaxy Spectra 
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[Fe/H] 
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[Co/Fe] 


[Ni/Fe] 
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1.94 


0.22 


6.2 


-0.07 


0.05 


0.02 


0.03 


-0.16 


0.05 


0.02 


0.03 


0.05 


0.03 


0.02 


0.02 


0.03 


-0.00 


2.05 


0.22 


6.5 


-0.05 


0.07 


0.08 


0.06 


-0.09 


0.08 


0.00 


0.03 


0.07 


0.01 


0.00 


0.03 


0.06 


-0.00 


2.14 


0.21 


6.8 


-0.04 


0.09 


0.13 


0.10 


0.02 


0.09 


0.03 


0.03 


0.08 


0.03 


-0.01 


0.03 


0.10 


0.01 


2.23 


0.22 


7.0 


-0.03 


0.13 


0.18 


0.17 


0.12 


0.12 


0.05 


0.03 


0.09 


0.03 


-0.01 


0.03 


0.16 


0.02 


2.30 


0.26 


7.6 


-0.01 


0.16 


0.21 


0.20 


0.22 


0.15 


0.09 


0.03 


0.11 


0.03 


-0.03 


0.02 


0.20 


0.03 


2.39 


0.29 


11.0 


-0.02 


0.19 


0.26 


0.25 


0.33 


0.20 


0.13 


0.04 


0.12 


-0.01 


-0.02 


0.05 


0.27 


0.01 


2.47 


0.34 


11.8 


0.00 


0.21 


0.27 


0.28 


0.43 


0.22 


0.16 


0.02 


0.12 


-0.02 


-0.03 


0.05 


0.26 


0.07 



Note. — The rms deviation between the model and data is computed over the 
less than 1% and so are omitted from the table. Systematic errors dominate the ei 

cent CMD-based de terminations (Feltzing & Johnson! 120021; 
Mom anv et al.M2003h . The elemental abundance data are de- 
rived from hi gh-resolution s pectra of individual RGB stars, 
as reported inlCarretta et all (12001b . i Zoc cali et al. (2004J), and 
OrigliaetaJJ (12005I) . 

The integrated light spectrum for M67 was constructed 
in ISchiavon et alj d2004[) by combinin g spectra of individual 
member stars with a Salpeter] (119551) IMF. M67 contains a 
large population of blue straggler stars, but these were omit- 
ted in the integrated spectrum. The spectrum covers the range 
3640A, to 5400A, at a resolution of 2.7 A. The S/N of the in- 
tegrated spectrum was not provide by Schiavon et al., so we 
adopt a nominal S/N of 100 per pixel. 

For NGC 6528 we use the integrated light spectrum 
from the library of globular cluster spectra presented in 
ISchiavon et alj (120051) . These spectra were obtained via drift 
scan observations, so the spectra reflect the true integrated 
light of the clusters. The spectrum covers the wavelength 
range 3360A to 6430A at a resolution of ~ 3.1 A. The S/N 
ranges from «50 per pixel at 4000 A to 150 per pixel at 
5000A. In this paper we analyze the spectrum extracted within 
the central 7". Several spectral regions had to be masked ow- 
ing to sub-optimal sky subtraction. These regions include 
4155A-4165A, 4535A-4565A, 4850A-4870A, and 5035A- 
5065A. 

For both spectra we use two wavelength intervals when 
fitting the model to these data: 4000A-4600A, and 4600A- 
5350A. 

4. TESTING THE MODEL AGAINST METAL-RICH STAR CLUSTERS 

Before discussing the main results of this paper, in this sec- 
tion we present the results of fitting our models to spectra of 
the metal-rich star clusters M67 and NGC 6528. Star clusters 
are of course simpler to model than galaxies, owing to the fact 
that they are at least approximately coeval and monometallic. 
They therefore can serve as useful tests of any SPS model. 

In Figure [7] we compare the best-fit ages, [Fe/H], and sev- 
eral a elements for these two clusters. Turning first to M67, 
we derive an age of 3 .5 ± 0.3 Gyr, in excellent agreement with 
the CMD-based age. We also find excellent agreement for 
Fe, Mg, Ca, and Ti. The abundances of O and Si are un- 
derestimated by 1 -2cr. Because the effects of O and Si on 
the spectrum are somewhat more subtle than the other a el- 
ements, higher quality integrated light spectra spanning the 
optical-NIR wavelength range for M67 are needed before a 
more definitive conclusion regarding these elements can be 
reached. 



11 wavelength range used in the fit. Ages are light-weighted. Form al statistical errors are 
: budget, and are estimated to be 0.05 dex or less. See Section l5~4l for further discussion. 

The model also provides an excellent fit to the age and 
[Fe/H] abundance of NGC 6528. The a elements also broadly 
agree with the literature data, although perhaps the most 
striking impression from Figure [7] is the large scatter be- 
tween the literature estimates for the abundance ratios. The 
most dramati c difference i s in th e estimated [Ca/Fe] abun- 
dance, where lOriglia et all (12005b find [Ca/Fe]= +0.4 while 
IZoccali et al] d2004 quote [ Ca/Fe]= -0-40 ; a dif ference of 
0.8 dex! We have plotted the lOriglia et alj (120051) results in 
grey bec ause these authors m odel high-resolution NIR spec- 
tra while lCarretta et all (120011) and lZoccali et al.l (12004) model 
high-resolution blue spectra. We also draw attention to the 
fact that these two authors find relatively modest enhance- 
ments in O and Mg, less than 0.2 dex. Thus, while NGC 6528 
is a-enhanced, it is not tremendously so. 

We do not compare results for C nor N because these ele- 
ments can be synthesized in the envelopes of red giant branch 
(RGB) stars, so abundances of these elements derived from 
such stars do not reflect the original (i.e., main sequence) 
abundances. Our model recovers the latter, so a direct com- 
parison with RGB-based abundances is not possible for C and 
N. Owing to the large discrepancy in reported literature abun- 
dances, even for important species such as Fe, O, Mg, Si, Ca, 
and Ti, we made no effort to compare to other elements as 
their derived abundances from high-resolution stellar spec- 
troscopy appear to be even more uncertain. 

Our model successfully recovers the ages, [Fe/H], and 
a element abundances of metal-rich star clusters. Further 
progress on testing SPS models must await 1) higher quality 
integrated light spectra of star clusters covering a wider wave- 
length range and at higher S/N than is presently available; 2) a 
thorough evaluation of the results from high-resolution stellar 
spectroscopy. Until a consensus emerges on the abundance 
patterns of the stars in metal-rich star clusters, such systems 
cannot be used to rigorously test abundance ratio predictions 
of SPS models. 

5. RESULTS 

We now turn to our main results — the analysis of high 
quality stacked spectra of early-type galaxies from the SDSS. 
In this section we present and discuss the overall quality of the 
fits, the detailed abundance patterns, and the derived effective 
temperatures of the stars. Table [T]presents the main derived 
parameters from modeling the stacked spectra. Results for 
the abundance ra tios of neutron cap ture elements Sr and Ba 
were presented in lConrov et alj (120 131) . and the results for the 
variation in the IMF with galaxy a will be presented in future 
work. We simply note here that the IMF trends with a are 
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Figure 5. Isochrones at 13.4 Gyr from the Dartmouth Stellar Evolution 
Database (DSEP) for a variety of chemical compositions. Top panel: 
Isochrones for a range in compositions approximately spanned by the early- 
type galaxies in our sample. The case with [Fe/H]= 0.0 and [a/Fe]= +0.4 is 
extreme, as no galaxies in the sample have such large a enhancements. No- 
tice that the isochrone with [Fe/H]= -0.1 and [a/Fe]= 0.2 is almost identical 
to the solar metallicity isochrone. Bottom panel: Same isochrones as in the 
top panel, with a constant shift in 7 c ff applied to the entire isochrone (the shift 
applied is shown in the legend). For the modest range in chemical composi- 
tions probed by our sample, the expected effective temperatures of the stars 
can therefore be well-approximated by a solar metallicity model with a single 
additional free parameter describing the offset in T^f. 

similar to those presented in lConrov & van D okkum (2012b). 

5.1. Overall Quality of the Fits 

We begin by discussing the overall quality of the model fits. 
In Figures [8] and [9] we show the data and best-fit model for 
the lowest and highest a bins. Within each panel the fluxes 
have been continuum-normalized by a high-order polynomial. 
These figures demonstrate the overall very high quality of the 
fits. In the lowest dispersion bin the rms between data and 
model is only 0.22%, and x 2 /dof= 1.1, indicating that our 
model is able to fully capture the observed variation in the 
spectrum to the precision allowed by the data. At the high- 
est dispersion bin, where the abundance ratios deviate most 
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Figure 6. Continuum-normalized stacked spectra of SDSS early-type galax- 
ies in three velocity dispersion bins. 

strongly from the solar-scaled values, the rms is 0.34%, with 
X 2 /dof= 4.7. It is not surprising that the fit is formally worse 
for the highest dispersion bin, both because the S/N is higher, 
which imposes more stringent demands on the model, and be- 
cause the strong deviations from solar-scaled abundance ra- 
tios means that the synthetic stellar spectra play an important 
role in the modeling. The synthetic spectra carry significant 
uncertainties (e.g., due to the treatment of the model atmo- 
spheres and incomplete line lists), so their greater weight in 
the model will tend to result in larger model uncertainties. 

In Figure [I0]we show the residuals now for all seven bins. 
The most striking feature of the residuals is their similarity 
across the sample. For example, the residuals around the 
Ca II triplet at 8500A-8700A are relatively significant, but the 
magnitude of the residuals doe s not vary strongly with galaxy 
velocity dispersion (see also Conroy & van Dokkum 2012b, 
where this was demonstrated on an object-by-object basis). 
This weak dispersion-dependence in the residuals is impor- 
tant for several reasons. First and foremost, it suggests that 
the trends in the derived parameters, to be discussed later in 
this section, are likely insensitive to remaining model system- 
atics. In addition, the fact that the residuals at high dispersion 
are similar to the lowest dispersion bin, where the abundance 
ratios are close to solar-scaled values, suggests that the model 
deficiencies lie in the base models constructed with empirical 
stellar spectra. It is worth pointing out that the MILES stellar 
library, which forms the core of the model in the optical, has 
typical S/N of ~ 150 A" 1 , so at the level of interest here the 
uncertainties in the empirical spectra probably play an impor- 
tant role (the residuals in Figure [l0]are smooth in part because 
the model and data have been smoothed to a = 350 km s -1 ). 

Some systematics with dispersion are evident, for example 
at 4500A and 4675A. These are probably due to features that 
are not well-captured in our synthetic spectral library. In fu- 
ture work we will investigate the sources of opacity in these 
wavelength regions in an attempt to further improve the mod- 
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Figure 7. Best-fit parameters for the star clusters M67 (top panel) and NGC 6528 (bottom panel). Results from our full-spectrum fitting of the integrated 
light (open blue circles) are compared to literature values based on high resolution spectra of individual member stars (red lines and symbols). Fiducial ages 
for M67 and NGC6528 ar e taken to be 4 and 1 1 Gyr, res pectively, based on isochrone fitting to their CMDs. Abundance patterns for NGC 6528 are from 
ICarretta et all EOOll COD. IZoccali et all <200l Z04) a nd Origlia eTall <2005l O05). For M67 we adopt a fiducial abundance pattern of 0.0 ± 0.1 (roughly la 
errors Tautvaisiene et al. 2000; Shetrone & Sandauist 2000). Overall the agreement between our model predictions and the literature values is good, although the 
large variation in the literature determination of many of the elemental abundances for NGC 6528 prohibits a more definitive test of our model. 



els. 



5.2. Abundance Patterns 



We now turn to the elemental abundance ratios derived from 
our sample of stacked SDSS early-type galaxy spectra. 

In FigureQT|we show the relations between velocity disper- 
sion and C and N (left panel), and the a elements O, Mg, Si, 
Ca, and Ti (right panel). In most cases the formal statistical 
errors are smaller than the symbol sizes. The abundance of 
both C and N increase from [C/Fe]«[N/Fe]w 0.0 at low a to 
0.2-0.3 at h igh a. This is in qualitative agreement with pre- 
vious results (iKelson et al . 2006; Graves et al. 2007; Schiavon 
l2007HSmith et al.ll2009h Uohansson et al.ll2012h . 

The derived abundance ratio patterns of the a elements 
show several interesting features. First, all ratios are close 
to the solar-scaled values at low a. Second, the trends with 
u are stronger for the lighter elements O and Mg compared 
to the heavier elements Si, Ca, and Ti. Mg appears to track 
O closely, with a slight preference for [Mg/0]< 0.0 at high 
a. Ca tracks Fe clo sely over the full sample, in agreement 
with previous work (ISaglia et al] 120021 ICenarro et al.l | 2003| 
Thom as et all l2003ah iGraves et all 120071; ISchiavonl 120071; 



Smith et all 120091; Uohansson et alj|2012l) . This suggests that 



a substantial amount of Ca is formed in Type la S Ne, as ex- 
pected from theoretical yields (No moto et al.lll984l) . We note 
that this is only the second reported measurement of Ti in 
early-type galaxies (following Johansso n et al.ll2012l) . and the 
second also of Si (following Worthev et al. 2 0131) . 

In Figure [12] we show the results for [Fe/H] in addition to 
the other iron peak elements V, Cr, Mn, Co, and Ni, all of 
which are measured here for the first time. The [Fe/H] abun- 



dance increases modestly over the sample by only 0.07 dex. 
This remarkable uniformity of the Fe abundance, in light of 
the substantial variation in the abundance of the light and a 
elements, poses a challenge to chemical evolution models of 
early-type galaxies. The abundances of the iron peak elements 
V, Cr, Mn, and Ni closely track Fe over the full sample, imply- 
ing that these elements form in the same nucleosynthetic sites. 
In stark contrast is Co, which increases in relative abundance 
to Fe as a increases, reaching [Co/Fe]« 0.27 at the highest 
dispersions. 

This result is explored further in Figure [13] where we show 
the abundance ratios of several elements with respect to O as 
a function of a. Clearly both N and Co track O remarkably 
closely. C also tracks O more closely than Fe, but [C/O] is not 
quite as constant as [N/O] and [Co/O], suggesting that C and 
N may not be forming in exactly the same environments. 

Although not shown in a figure, the [Na/Fe] abundance also 
varies considerably across the sample, from -0.16 to 0.43, 
qualitatively tracking the light a elements (see Table [TJ. In 
detail however the Na abundance variation is considerably 
stronger than the other elements. The only two Na lines avail- 
able to us are the Na I features at 5895 A and 8190A, re- 
spectively. The former is well-known to be affected by the 
presence of an interstellar medium (ISM). However, the sam- 
ple was selected to have no detectable emission in Ha nor 
in [O II], suggesting that contamination by ISM absorption 
should be small. Moreover, one might have expected ISM 
contamination to be greater at lower cr, because low mass 
galaxies are younger, and yet the derived [Na/Fe] at low a is 
lower than either [O/Fe] or [Mg/Fe] (the presence of an ISM 
would cause stronger absorption at 5895A and thus larger 
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Figure 8. Comparison between the SDSS early-type galaxy stacked spectrum in the lowest velocity dispersion bin (black lines) and the best-fit model (red lines). 
Within each wavelength range, the top panel shows the continuum-normalized fluxes and the lower panel shows the ratio between data and model. The grey band 
denotes the Icr error limits of the data. Overall the quality of the fit is excellent. 

derived [Na/Fe] values). Of course, it is not necessary that 
[Na/Fe] track the other elements in detail, as its nucleosyn- 
thetic origins are quite complex, including contributions from 
both massive stars and intermediate-mass asymptotic giant 
br anch (AGB) stars. In deed, the chemical evolution model 
of lArrigoni et al.l (120101) finds a variation in [Na/Fe] in broad 
agreement with what we find herein, ranging from -0.2 at low 
mass to +0.3 at high mass. 



5.3. Effective Temperatures 

As discussed in Section [2] one of the free parameters in 
our model is the shift in T e s from a fiducial solar metallicity 
isochrone. This is an important parameter to include because, 
as we saw in the previous section, the metallicities and abun- 
dance patterns vary by factors of several across our sample. 
A change in metallicity and/o r abundance patte rn will induce 
a change in T e ff of the stars (Pott er et al.ll2007h . It would be 
difficult to include this variation self-consistently owing to the 
large grid of isochrones that would be required. The parame- 
ter Ar e ff is meant to capture all of these effects into a single 
number. 

In Figure [14] we show the resulting relation between AT e ff 
and the sum of [O/Fe] and [Fe/H], which is a proxy for [Z/H] 



(e.g.. lTrager et al.ll2000b . The relation derived for our sample 
of SDSS early-type galaxies is compared to theoretical expec - 
tations from the Dartmouth isochrones dDotter et al.l 2008). 
The early-type galaxy sample increases in a from low to high 
metallicity. The theoretical predictions are based on the mod- 
els shown in Figure [5] The star symbol corresponds to the 
model with [Fe/H]=-0.1 and [a/Fe]=+0.2 and highlights the 
fact that Ar e ff is not a simple function of [Z/H]. 

The excellent agreement in the slope of the AT e s vs. 
[0/Fe]+[Fe/H] relation between the theoretical prediction and 
the data suggests that we are reliably measuring the distribu- 
tion of effective temperatures for the stars in these galaxies 
(i.e., we are solving for the detailed shape of the isochrones). 
The flattening of the observed trend at high a is interesting 
and will be explored in detail in future work. The offset be- 
tween the two relations is ps 40K, which is well within both 
the systematic uncertainties in the color-temperature relations 
used to assign temperatures to empirical stars and the theoret- 
ical uncertainties in modeling RGB temperatures. 

In order to provide a more direct comparison between the 
observations and the theoretical predictions we will need to 
create isochrones with the exact same chemical composition 
as we derive for the early-type galaxies. This will be the sub- 
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Figure 9. Same as Figure[8] now showing results for the highest velocity dispersion bin. As in Figure|8]the overall quality of the fit is excellent. 



ject of future work. 

5.4. Exploration of Systematic Uncertainties 

In this section we explore variations in the fiducial model 
and in our primary set of stacked spectra in order to quan- 
tify systematic uncertainties in our main results. In Figure [TBI 
we plot the difference in either abundance ratio or log(age) 
between our main results and a modification to our standard 
analysis. In the upper left panel we compare our standard 
analysis to a simplified model with fewer free parameters. In 
particular, in this model the IMF is fixed to the Kroupa (2001 ) 
form, only a single age component is included (rather than the 
standard two components), the iron peak elements are forced 
to track Fe, and all the nuisance parameters are turned off (in- 
cluding Ar e ff). This model results in a larger rms residual 
between data and model compared to our fiducial model (i.e., 
it is a poorer fit), which should not be surprising given that 
there are fewer free parameters. Nonetheless, the resulting 
abundance ratios agree well with the standard model, in par- 
ticular Mg, C, Ca, Si, Ti, and Fe agree to within 0. 1 dex and 
in most cases to within 0.05 dex. The largest outliers are the 
age at low a and the [O/Fe] and [N/Fe] abundances. The dif- 
ference in [O/Fe] is particularly strong. We speculate that the 
larger range in [O/Fe] in the simple model is due to the fact 
that the parameter Ar e ff is not included. In essence, [O/Fe] 



is forced to vary more in order to compensate for the lack of 
variation in AT e ff. 

In the upper right panel we consider a model where the 
mean wavelength-dependent residuals averaged over the sam- 
ple are added back into the model. In other words, in this case 
the model is multiplied by the mean residuals shown in Figure 
[l0]before the model is fit to the data. This will obviously re- 
sult in a much lower rms difference between the best-fit model 
and data, and the question is whether or not this changes the 
derived parameters. As demonstrated in Figure [15] the re- 
sulting change in abundance patterns is very small, typically 

< 0.03 dex. This implies that the small mis-match between 
the data and our model does not bias the resulting best-fit pa- 
rameters. 

In the lower left panel we consider the change in parameters 
when only the blue spectral region (A < 5800A) is included 
in the fit. Here again most parameters change very little, by 

< 0.05 dex. Significant outliers include [Na/Fe] and [O/Fe]. 
The change in [Na/Fe] is not surprising because the two fea- 
tures most sensitive to Na are the Na I features at 5895A and 
8190A. With these features masked out the constraint on the 
Na abundance is coming exclusively from the effect of Na 
on the electron pressure in the stellar atmospheres, which in- 
directly affects the ionization states of other elements. For 
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Figure 10. Residuals between the data and best-fit models for all seven stacked early-type galaxy spectra analyzed in this paper. The shaded region marks the 
noise limit of the middle <r bin. In most spectral regions the residuals are constant across the sample, which is an important result because it means that the 
higher a galaxies, which harbor abundance patterns significantly different from the solar neighborhood, are nearly as well described by our model as the lower 
a, less metal-enhanced galaxies. In a few spectral regions, e.g., at 0.45/im the residuals correlate with galaxy velocity dispersion, indicating that our model is not 
capturing the totality of variation in the observed spectra. 



[O/Fe] the deviation is also significant, exceeding 0. 1 dex in 
several bins. The origin of this offset is not clear, although we 
note that the best-fit AT^g parameter is lower in the case of 
fitting only the blue spectra. With cooler stars the metal-lines 
are stronger, perhaps alleviating the need for higher [O/Fe] 
values. This highlights the need for spectra that extend into 
the red, where the Ar e ff parameter has a large effect and hence 
can be more robustly measured. Moreover, at A < 5800A the 
constraint on [O/Fe] is coming primarily through the effect of 
O on molecular dissociation equilibrium in the stellar atmo- 
spheres, while at A > 5800A the effect of O on the spectrum 
is primarily to increase the strength of the TiO lines (see Fig- 
ures [2] and |3j. The red spectrum therefore offers at the very 
least a complementary, and arguably a stronger constraint on 
the [O/Fe] abundance. 

In the lower right panel we show the results of fitting to 
the unsmoothed stacked spectra. This is primarily a test of the 
model in the low a bins, where the unsmoothed spectra have a 
resolution w4x higher than the smoothed spectra. Here again 
essentially all parameters agree with our standard results to 
better than 0.03 dex. The largest outlier at low a is [O/Fe], 
but even in this case the disagreement is < 0. 1 dex. 

We have performed an additional test (not shown in the fig- 
ure), where we masked wavelengths bluer than 4200A. This 
was done primarily to test the [Co/Fe] result, since there are 
several strong Co I lines with hyperfine splitting components 



at A < 4200 A and many weaker blends at longer wavelengths 
(see Figure |2). The differences in the derived [Co/Fe] ratios 
between this and our fiducial model is < 0.03 dex for all a 
bins except for the lowest bin, where the difference is 0.046 
dex. The result for Co seems be insensitive to the particular 
wavelength range used in the fit. 

In summary, the tests performed in this section suggest that 
the systematic uncertainties in our derived parameters is prob- 
ably < 0.05 dex. 

6. COMPARISON TO PREVIOUS WORK 

In this section we provide a comparison to other tech- 
niques for deriving stellar population parameters of early - 
type galaxies. In Figure [16] we show abundance ratios and 
light-weighted ages from several groups, including the results 
presented here, from the EZ Ages code (Graves & Schiavon 

•lotr 

W13). 

The EZ _Aqes code fits 7 Lick indices with the SPS mod 



and fro m lThomas et al.1 (201 
. J12), and lWorthev eTafl (120131 



[HI T10), Uohansson et al 

m 



els of ISchiavonl d2007h in order to derive the age, [Fe/H], 
[C/Fe], [N/Fe], [Mg/Fe], and [Ca/Fe]. The code was run on 
the exact same stacked spectra analyzed herein, so that any 
differences can be attributed entirely to differences in model- 
ing techniques. In general the agreement is very good. The 
[C/Fe] and [N/Fe] trends are stronger with EZ_Ages com- 
pared to the results presented heren. Perhaps the most signifi- 
cant differences are the overall higher [Mg/Fe] abundances in 
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Figure 11. Abundance ratios of C, N, and the a elements as a function of early-type galaxy velocity dispersion. In most cases the statistical errors are smaller 
than the symbol sizes. C, N, O, and Mg all seem to track each other fairly closely, w hile the heavier a elements Si, Ca, and Ti show weaker trends with a. Ca 
clearly tracks Fe instead of O. Systematic errors are probably < 0.05 dex (see Section [5~4V 




Figure 12. Abundance ratios of the iron peak elements as a function of early-type galaxy velocity dispersion. In most cases the statistical errors are smaller than 
the symbol sizes . All iron peak elements track Fe except for Co, which seems to track the light a elements O and Mg. Systematic errors are probably < 0.05 dex 
(see Section l5~4l . 



EZ_Ages (by w 0.07 dex), and the lower [Fe/H] abundances 
(by « 0.07 dex). 

The results from T IP a re based on t h e SPS models of 
iThomas et all (l2003bh and iThomas et all (120041) . T10 ana- 
lyzed a morphologically-selected sample of early-type galax- 
ies from SDSS. They used 25 Lick indices in order to de- 
rive three parameters: [Z/H], age, and [a/Fe]. We compare 
their derived [a/Fe] values to both [Mg/Fe] and [O/Fe]. We 
find overall satisfactory agreement between their [a/Fe] val- 
ues and our [Mg/Fe] and [O/Fe] values, and also between their 
ages and ours. 

J 12 present results from the exact same morphologically- 
selected sample of early-type galaxie s analyzed in TIP, bu t 
with the updated SPS models from Thomas et al] d2Pl lb . 
These models allow for variation in the elements C, N, O, 
Mg, Ca, and Ti, in addition to metallicity and age. They fit 



their models to the data using 18 Lick indices. The differ- 
ences between J 12 and our own work is the largest of any of 
the techniques shown in Figure [16] (except for [O/Fe] from 
W13), although the differences rarely exceed P.l dex. It is 
difficult to know if the differences are due to modeling tech- 
niques or sample selection, especially at low a where there is 
likely to be significant scatter between morphological ('early- 
type') and emission line ('quiescent') selections. 

Very recently, W13 analyzed stacked SDSS early-type 
galaxy spectra using their SPS models. The stacks are not 
identical to those analyzed herein, but are very similar and so 
a direct comparison between W13 and our results is possible. 
From Figure Q~6]it is clear that the ages, [Mg/Fe], [N/Fe], and 
[Ca/Fe] trends are in good agreement between the techniques, 
except for the age in the highest a bin in W13. The [C/Fe] 
abundances are also in reasonable agreement. In contrast, the 
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Figure 13. Abundance ratios between C, N, Co, and O as a function of early- 
type galaxy velocity dispersion. All of these elements seem to track each 
other very closely, espe cially N, O, and Co. Systematic errors are probably 
< 0.05 dex (see Section rllt . 
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Figure 14. Relation between offset in isochrone T e ff and metal abundance. 
Results are shown for the stacked SDSS early-type galaxy data and com- 
pared to expectations based on theoretical isochrones from the Dartmouth 
Stellar Evolution Database (DSEP). The theoretical offsets are derived from 
Figure [5] The star symbol corresponds to the model with [Fe/H]= —0.1 and 
[a/Fe]=+0.1 and highlights the fact that Ar c ff is not a simple function of 
[Z/H]. The dotted lines are offset from each other by 40K. The excellent 
agreement in the slope between our modeling and the DSEP predictions is a 
strong indication that our model is able to reliably measure the distribution of 
r c ff in early-type galaxies. 

[O/Fe] trend derived by W13 is much steeper than not only 
our results but also those of J12 (and T10 if one treats their 
[a/Fe] as [O/Fe]). In addition, their inferred [Fe/H] abun- 
dance^ range from -0.05 at low a to -0.7 at high a. The 

6 W 13 (in addition to T10 and J 12) do not constrain [Fe/H] directly but in- 
stead estimate [Z/H] and the abundance ratios. [Fe/H] is thus a derived prod- 
uct and can be roughly estimated via [Fe/H]=[Z/H]-A [a/Fe] with A in the 
range 0.7 — 0.95 depending on the detailed abundance pattern fTrager et all 
2000). For the purposes of this discussion we estimated [Fe/H] from W13's 
data via [Fe/H]=[Z/H]-0.9[O/Fe]. 
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Figure 15. Effect on derived parameters due to varying one or more aspects 
of our analysis (plotted here as fiducial model minus modified model). All pa- 
rameters in TablefT]are plotted here (except in the top left panel), but only the 
most deviant parameters are highlighted. Top Left Panel: Variation between 
the standard model and a simple model that contains many fewer parameters. 
Top Right Panel: Variation between the standard model and one in which 
the mean residuals (averaged over all <r bins) are added back into the model. 
Bottom Left Panel: Variation in parameters when fitting only the blue spectral 
region (A < 5800A). Bottom Right Panel: Variation in parameters when the 
spectra are not smoothed to a common dispersion of <r = 350 km s _1 . The sin- 
gle error bar in each panel represents an error of ±0.05 dex, a value we regard 
as a plausible upper limit to the systematic errors on the derived parameters. 

results for [Fe/H] are so different from the other results in 
Figure [16] that they were omitted from the Figure for clarity. 
W13 fixed [Ti/Fe]= 0.0 throughout their analysis. These au- 
thors also presented results for Si, finding [Si/Fe] is approxi- 
mately 0.0 at low a and increases to 0.5 at high a. This is also 
in contrast to our results, in which [Si/Fe] increases from 0.0 
to 0.16 from low to high a. At present we do not know the 
origins of these discrepancies. 

One of the more notable features of Figure [16] is the 
[Mg/Fe] vs. a panel, where there is a ^0.1 dex variation 
in the normalization of the relation. This variation is signif- 
icant in the context of using [Mg/Fe] t o infer a star forma- 
tion timescale (e.g.. [Thomas et alfcooBI) . iTrager et all (120081) 
noted that their models produced [Mg/Fe] values lower than 
previous Lick index-based models, and they speculated that 
this may be due to the different treatment of the response 
functions. In all previous work on Lick index modeling, three 
theoretical stellar spectra were combined to estimate the re- 
sponse of the integrated ligh t spectrum to abu ndance vari- 
ations dTripicco & Belli 119951: iKorn et all 120051) . Trager et 
al., and also Worthe y et al.1 (120131) . instead used many stellar 
spectra covering the isochrone in order to estimate response 
functions, in a manner similar to what is done in the CvD 
models. In addition, previous Lick index-based models cal- 
culated response functions for an isochrone with an age of 5 
Gyr, whereas more recent models use older isochrones (or a 
range of isochrones appropriate for the ages of the models). 

We have estimated the effects of these assumptions by mod- 
ifying our model to employ only 3 stars (as opposed to the 
usual 20) to estimate response functions, and also to use stars 
along a 5 Gyr rather than a 13 Gyr isochrone. The number 
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Figure 16. Comparison between different t echniq ues for estim ating ages and eleme ntal ab undances. Our results are compare d to techniques that rely on Lick 
index fitting, including lGraveT & Schiavon (2008. EZ_Ages), Thomas eTaTl <2010l T10), Joha nsson et"aT] <20T1 J12), and lWorthev"eTaI] i2013l W13). The 
EZ_Ages results are based on modeling the exact same stacked spectra as used herein, the W13 results are based on similar but not identical stacks, while the 
T10 and J12 results are linear fits to results obtained for ~ 4000 SDSS early-type galaxies with light-weighted ages greater than 2.5 Gyr. Note that T10 derive 
a single [a/Fe] ratio, which we plot in both the [Mg/Fe] and [O/Fe] panels, as it is not entirely clear which element their [a/Fe] parameter is most sensitive to. 
Qualitatively the derived trends agree well between various methods/groups, although quantitative differences are clearly evident. 



of stars used to compute response functions has a significant 
effect. We have created response functions for the Mg b Lick 
index due to a factor of two change in [a/H]. A model using 
20 stars has a response function in the Mg b equivalent width 
(EW) that is 30% greater than a model that uses only 3 stars. 
The effect of the age is rather modest, with a change in EW of 
only 5%, with the 13 Gyr isochrone producing greater EWs. 
The combination of both effects is a 37% increase in the re- 
sponse function of the Mg b index when 20 stars are used 
along a 13 Gyr isochrone, compared to only 3 stars along a 
5 Gyr isochrone. This translates directly into differences in 
the derived [Mg/H] abundances, and therefore any model us- 
ing only a small number of stars along an intermediate-age 
isochrone will tend to overestimate the [Mg/H] abundance. 

In summary, there is broad agreement between models for 
those elements that impart strong, localized changes in low 
resolution data (especially C, N, Mg, and Ca; see Figure [2). 
In such cases models based on spectral indices seem to per- 
form as well as our model based on full spectrum fitting. The 
power of the full spectrum fitting technique lies in measur- 



ing parameters that either impart a more subtle change in the 
spectrum (as in the case of Si and Mn), or affect large swaths 
of the spectrum (as in the case of O, Ti, and most of the iron 
peak elements). 

7. DISCUSSION 

The main results from this paper are summarized in Fig- 
ure[17] where we show the abundance patterns of 13 elements 
as a function of early-type galaxy velocity dispersion. This 
is the most comprehensive analysis of the detailed abundance 
patterns of the stars within distant galaxies to-date. We em- 
phasize that the spectra sample approximately the inner Q.5R e , 
and thus the abundances derived should be representative of 
the inner regions of the galaxies. In this section we discuss 
several implications of these results. 

The qualitative agreement between our results and those of 
previous authors is, for the most part, very encouraging. This 
is significant because our model is a major departure from 
all previous analyses of the abundance patterns of early-type 
galaxies. Previous work focused on the analysis of selected 
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Figure 17. Summary of the abundance trends derived herein from stacked spectra of SDSS early-type galaxies. The sample spans a range in velocity dispersion 
from cr = 90 kms~' to a = 300 kms~' . The results for Na are omitted from this summary diagram, though we note that [Na/Fe] varies from —0.16 to +0.43 from 
low to high cr. 

ratios that we derive herein have important implications for 
the inferred formation timescales of the stars within these 
massive galaxies. As an example of the effect, Thom aset alj 
( 120051) use the equation [a/Fe]« | - i l og At, calibrated from 
the simple chemical evolution model of Tho mas et al.l ( 1999), 
in order to convert their measured [a/Fe] ratios into con- 
straints on the formation histories of early-type galaxies. In 
our analysis we find [Mg/Fe]= 0.22 in the most massive galax- 
ies, which, when employing the equation above, corresponds 
to a timescale of 0.8 Gyr. In contrast, the EZ_Ages model 
and T10 find [Mg/Fe]w 0.27 at the highest dispersions, corre- 
sponding to a timesc ale of 0.4 Gyr. Based on the models of 
IThomas et all (120111) , J12 find [Mg/Fe] « 0.33 at the highest 
dispersions, corresponding to a remarkably short timescale of 
0.2 Gyr. It is interesting in this context that the [O/Fe] ra- 
tios derived by J 12 are in much better agreement with our 
results at the highest dispersions, and so if [O/Fe] were used 
as a tracer of a, rather than Mg, then the inferred timescales 
would be more similar between different groups at high a. 
However, this of course implies that [O/MgJ^O, and in this 
case it is not clear that a single parameter "[a/Fe]" can be 
used to simply (or reliably) translate derived abundance ratios 
into SF timescales. Clearly sophisticated chemical evolution 
models are need to provide further insight. 

One of the more interesting trends evident in Figure [17] 
(see also Figure [T3l is the fact that the iron peak element Co 
closely tracks O rather than Fe. This is the first measurement 
of Co in the integrated light spectra of early-type galaxies, and 
it provides an important constraint on the nucleosynthetic ori- 
gin of Co. Taken at face value, [Co/Fe]«[0/Fe] implies that 
Co forms primarily in massive stars, in apparent contrast to 
the other iron peak elements. In this context it is worth recall- 
ing that Co is also peculiar am ongst the iron pea k elements 
in the metal-poor Galactic halo (Mc William 1997). This may 
be a coincidence at some level, as the abundance pattern of 
the very metal-poor halo is probably due exclusively to very 
metal-poor Type II SNe, in contrast to the processes operating 



spectral indices with theoretical spectral models that were rel- 
atively simple with respect to the theoretical spectral libraries 
used herein. In addition, we have included a large number of 
additional free parameters in order to marginalize over our ig- 
norance of various aspects of stellar evolution and stellar pop- 
ulations. In spite of these differences in techniques, there is 
excellent agrement in the trends, including the derived light- 
weighted ages, the weak variation in [Fe/H] and [Ca/Fe], and 
the strong variation in [Mg/Fe], [C/Fe], and [N/Fe]. It should 
not be surprising that there is good agreement between tech- 
niques for the elements C, N, Mg, and Ca because these ele- 
ments impart strong, localized changes in low resolution spec- 
tra (see Figure |2}. Quantitatively there are differences in the 
slope and normalization for some of these relations, though 
the differences rarely exceed 0.1 dex in amplitude. Based on 
our preliminary investigations, it appears that at least some of 
the differences can be traced back to the number of theoretical 
spectra used to estimate re sponse functions to elemental abun- 
dance variations (see also Tra ger et al.l l2008). The real power 
of the full spectrum fitting technique is in the measurement of 
parameters that either impart a small change in the spectrum 
(e.g., Mn, Sr, Ba), or that affect large regions of the spectrum 
(e.g., O, Ti, V, Cr, Mn, Fe, Co, Ni). 

The standard interpretation of [a/Fe] is that this ratio is sen- 
sitive to the time-scale of star formation (SF), with higher val- 
ues corresponding t o shorter timescales (e.g., Tinslevl ll979l : 
IThomas et alj [l999). As mentioned in the Introduction, this 
ratio is also sensitive to the IMF, the nucleosynthetic yields, 
the delay time distribution of Type la SNe, and the preferen- 
tial loss of metals via winds. Some chemical evolution mod- 
els are able to reproduce the observed variation in not only 
Mg and O but also Si and Ca based on the assumption that the 
main variable drivin g these trends is indeed the SF timescale 
(Pipin o~et al.l 12009), while others also require variatio n in 
other parameters such as the IMF ( Arrigoni et al.lEOlOl) . In 
the context of any chemical evolution model, the [Mg/Fe] ra- 
tio will be correlated with the SF timescale, and so the lower 
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to shape the abundance patterns in massive early-type galax- 
ies. Core-collapse SNe nucleosynthesis models predict that 
Co is synthesized by complete Si-burning in the deepest lay- 
ers, while Cr and Mn are produced in the outer incomplete 
Si-bu rning layers (Woosley & Weaver 1995; Nak amura et al.l 
1999). There is thus a mechanism by which to separate the 
production of Co from the other iron peak ele ments. Perhaps 
more s ignificantly, the Type la SNe yields from lNomoto et all 
( 1984) show a deficit of Co production compared to other iron 
peak elements. Again, detailed chemical evolution models 
will be required to sort out these details. 

We reaffirm previous results that the light a elements (e.g., 
O and Mg) show greater variation than the heavier a elements 
(e.g., Ca and Ti). We also find in this work that the variation 
in Si is intermediate between the lighter and heavier elements. 
At least some of these tren ds can be understood i n the context 
of Type la SNe yields. The Nomoto et al] (119841) models pro- 
duce a large quantity of Ca (comparable to Cr, Mn, Fe, and 
Ni), and much more Si than either O or Mg. The modest vari- 
ation in Ti is somewhat more difficult to understand in this 
context, because the Nomoto et al. models produce very little 
Ti (comparable to O and Mg). Whether or not the variation in 
[Ti/Fe] with galaxy dispersion can be reproduced with stan- 
dard nucleosynthesis yield tables remains to be seen. 

Finally we turn to C and N. As found in previous work, 
both [ C/Fe] and [N / Fe] increase with increasing disper- 
sion dSchiavonl|2007HGraves et al.ll2007t ISmith et all 120091; 
iJohansson et al. 120121) . We find that both of these elements 
track O and Mg quite closely (in detail C seems to track Mg 
while N tra cks O, but we regard this level of similarity as for- 
tuitous). Johans son et al.l(l2012l) also find that C closely tracks 
O, with [C/0]~ 0.0, but they find sub-solar [N/O] values. The 
fact that C and N are tracking the light a elements is probably 
telling us either that C and N form mostly in massive stars, or 
that there is a lower limit to the star formation timescale such 
that both massive stars and moderately massive AGB stars 
(which can produce lots of C and N) are always contributin g 
to the enrichment of the ISM (see also Johan sson et al.ll2012l) . 
The latter interpretation is consistent with the star formation 
timescales derived using our [Mg/Fe] ratios. Much effort is 
still needed on the theoretical side to understand the nucle- 
osynthetic origins of these elements in order to use them as 
cosmological clocks. 

8. SUMMARY 

In this paper we have used the CvD population synthesis 
model to fit very high quality stacked spectra of SDSS early- 
type galaxies spanning a range in velocity dispersion from 
90 kms -1 to 300 kms -1 . The novel feature of our approach 
is fitting the full continuum-normalized optical-NIR spectra 
(4000A-8800A), rather than a select number of indices. The 
model includes variation in age, the abundances of 16 ele- 
ments, the effective temperature distribution of the stars, the 
stellar IMF, amongst other parameters. The resulting qual- 
ity of the fits are very high, with rms residuals ranging from 
0.2% -0.3% across the sample. A variety of tests reveal that 
the systematic uncertainties in our measurements are probably 
0.05 dex or less. We now summarize our main results. 

• Across the sample from low to high a, ages increase 
from 6 to 12 Gyr, [Fe/H] increases from -0.07 to 0.00, 
and the abundance ratios [O/Fe], [Na/Fe], [Mg/Fe], 
[Si/Fe], [Ti/Fe], and [Co/Fe] increase from approxi- 
mately zero to 0.1 -0.4, with O, Na, Mg reaching the 



highest values. In contrast, [Ca/Fe], [V/Fe], [Cr/Fe], 
[Mn/Fe], and [Ni/Fe] remain nearly constant as a in- 
creases. 

• N and Co track O to better than 0.03 dex over the full 
sample, and C tracks O to better than 0. 1 dex. Co there- 
fore probably forms predominantly in massive stars. 
For N and C the situation is more complex because 
these elements probably form in both massive stars and 
AGB stars. The constant values of [C/O] and [N/O] 
therefore either requires C and N to form mostly in 
massive stars or the star formation timescale to never 
be shorter than several hundred Myr, the lifetime of a 
~ 3M Q star. 

• Our model allows us to measure the shift in the effective 
temperature of the stars with respect to a solar metallic- 
ity isochrone. This shift inversely correlates with the 
total metallicity, with a slope that is in excellent agree- 
ment with theoretical expectations. We are therefore 
able to measure not only the mean ages and abundance 
patterns but also the distribution of effective tempera- 
tures of the stars in early-type galaxies. 

• Comparison of our results to Lick index-based tech- 
niques reveals agreement in several quantities, such as 
age, [N/Fe], and [Ca/Fe], and less agreement in others 
including [C/Fe], [Fe/H], and [Mg/Fe]. Some of the 
disagreement may be due to different samples of early- 
type galaxies, but the different modeling techniques 
also seems to play an important role. Perhaps most 
relevant from the standpoint of broader implications, 
we find lower [Mg/Fe] values at fixed a (by 0.05-0.1 
dex) than s ome previous wor k (especially those which 
employ the iKorn et alj 12005, response tables), which 
we believe is due to our more accurate modeling of 
spectral variations due to abundance changes. This 
implies longer formation timescales for massive early- 
type galaxies than previously reported. 

These results will set the foundation for future studies 
aimed at constraining the evolutionary histories of galaxies 
from their detailed abundance patterns. With such a large 
number of elements precisely measured over a wide range in 
galaxy mass, it should be possible to place novel constraints 
not only on the formation histories of these galaxies but also 
on the nucleosynthetic origins of elements that are difficult to 
model from first principles. 
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